Soliton ratchets in homogeneous nonhnear Klein-Gordon systems 



Luis Morales-Molina 

Max-Planck Institut fiir Physik Komplexer Systeme, Nothnitzer Sir. 38, 01187 Dresden, Gejman^ 

Niurka R. Quintero 

Departamento de Fisica Aplicada I, E.U.P., Universidad de Sevilla, c/ Virgen de Africa 7, 41011 Sevilla, Spaing 

Angel Sanche^ 

Grupo Inter disciplinar de Sistemas Complejos (GISC) and Departamento de Matemdticas, 
Universidad Carlos III de Madrid, Avenida de la Universidad 30, 28911 Leganes, Madrid, Spain, and 
Institute de Biocomputacion y Fisica de Sistemas Complejos (BIFI), 
Universidad de Zaragoza, 50009 Zaragoza, Spam 

Franz G. Mertenfi 
Physikalisches Institut, Universitdt Bayreuth, D-85440 Bayreuth, Germany 
(Dated: February 2, 2008) 

We study in detail the ratchet-like dynamics of topological solitons in homogeneous nonlinear 
Klein-Gordon systems driven by a bi-harmonic force. By using a collective coordinate approach 
with two degrees of freedom, namely the center of the soliton, X{t), and its width, l{t), we show, 
first, that energy is inhomogeneously pumped into the system, generating as result a directed motion; 
and, second, that the breaking of the time shift symmetry gives rise to a resonance mechanism that 
takes place whenever the width l{t) oscillates with at least one frequency of the external ac force. 
In addition, we show that for the appearance of soliton ratchets, it is also neccesary to break the 
time-reversal symmetry. We analyze in detail the effects of dissipation in the system, calculating 
the average velocity of the soliton as a function of the ac force and the damping. We find current 
reversal phenomena depending on the parameter choice and discuss the important role played by the 
phases of the ac force. Our analytical calculations are confirmed by numerical simulations of the full 
partial differential equations of the sine-Gordon and 4>^ systems, which are seen to exhibit the same 
qualitative behavior. Our results show features similar to those obtained in recent experimental 
work on dissipation induced symmetry breaking. 

LEAD PARAGRAPH 

Ratchet, or rectification, phenomena in nonlinear nonequilibrium systems are receiving a great deal of attention in 
the past few years. A typical example of a ratchet occurs when point-like particles are driven by deterministic or 
non-white stochastic forces: Under certain conditions related to the breaking of symmetries, unidirectional motion can 
take place although the applied force has zero average in time. The broken symmetries can be spatial, by introducing 
an asymmetric potential, or temporal, by using asymmetric periodic forces. This effect has found many applications 
in the design of devices for rectification or separation of different particles. Recently, ratchet dynamics is being studied 
in the context of solitons, trying to understand whether the fact that solitons, or nonlinear coherent excitations in 
general, are extended objects, allows for similar rectification phenomena. In nonlinear Klein-Gordon systems, it has 
been shown that ratchet-like behavior can be observed when driven by asymmetric periodic forces. However, the 
phenomenon is different from that seen in point-like particles as the deformations of the soliton play a crucial role. 
In this paper we describe in detail these phenomena, analyze the mathematical conditions for its existence, explain 
the underlying mechanism originating the net motion, and study the highly non-trivial effects of dissipation on the 
motion. Our results may be related to recent observations of rectification in Josephson junctions and optical lattices. 

I. INTRODUCTION 

Ratchet or rectification phenomena, the subject of an intense research effort during the last decade, appear in very 
many different fields ranging from nanodevices to molecular biology 
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appearance of ratchet-like behavior requires two ingredients: Departure from thermal equilibrium, either by u sing 
correlated stochastic forces or deterministic forces, and breaking of spatial and/or temporal symmetries pL la. IH flOl 
im IT^ IT^ IT3 | . Many applications of the ratchet theoretical framework deal with the scenario of spatial symmetry 
breaking, beginning with the proposal for nonequilibrium rectifiers in T|. However, in other problems, breaking the 
temporal symmetry may be much more feasible; this is the case, for instance, of optical lattices [l5l | or Josephson 
j mictions m llil (see also ^ for related work on Josephson ratchets in the dissipative quantum regime) when 
driven by an appropriate bi-harmonic force. In this context, the ratio of the frequencies of a bi-harmonic force is 
related with the breaking of the temporal-shift symmetry, ( give n by f{t) = —f{t + T/2) with T being t he p eriod) 
whereas suitable choices for the driving phases and damping lead to time-reversal symmetry breaking [2l|. This 
mechanism to obtain directed motion from a zero-mean force has been first proposed for particle (zero-dimensional) 
systems 0,11,01' ^^'^ ^^^.s recently been extended to extended systems, both classical 0,0,123, HI] and quantum 
ones Q|. 

A very relevant system, mostly because its character of paradigmatic of phenomena arising in connection with 
topological solitons and its many applications, is the damped and driven sine-Gordon (sG) model |25i |: 

- (l>xx + sm{(j>) = -(3(l)t + fit). (1) 

In the symmetry considerations above were studied in the context of Eq. ^ by choosing 

f(t) = ei sin{5t + So) + €2 sm{mSt + So + 9') = ei sin{St + Oi) + £2 sin(TO(5t + 6*2), (2) 

where m is an integer number and ei (£2) is the amplitude of the harmonic component with frequency S (rnS) and 
phase 9i {62) ■ We note that in the equation above 9i = So y d' = 92 — 9i, being 9' the relative phase, and therefore the 
results presented here can be understood as well in terms of relative phases. The reason for including two harmonics 
is that, as can be seen immediately from the previous symmetry analysis, in the case of a single harmonic the directed 
motion of sG kinks is not possible, something advanced by a different reasoning some years ago |26j and experimentally 
confirmed in [23|. Thus, the main result in 2.3j was the finding that, if f{t) breaks the shift symmetry and if the 
total topological charge in the system is nonzero^ a directed current (meaning a preferential motion of the kink in 
one direction) could be observed whose direction and magnitude depend on the driving parameters and damping 
coefficient. 

Subsequently, in a recent paper l2^ , we advanced a first explanation of the mechanisms underlying the observations 
in pS] : Directed motion of solitons arises because of a resonance phenomenon involving the oscillations of the soliton 
width, l(t), with at least one of the frequencies of the ac force. In terms of the symmetry picture we are discussing 
here, we note that it is precisely the breaking of the shift symmetry the reason for the appearance of the resonance, 
which means that the shift symmetry and the resonance conditions are equivalent. In the present work we will extend 
our preliminary results in 28] in several directions. First, we will show that net motion of the soliton is only possible 
if, in addition to the aforementioned resonance condition on the frequencies, the time-reversal symmetry is also broken 
by the action of ac force and damping. As will be seen below, this in turn implies that the energy introduced by 
the biharmonic force is inhomogeneously pumped into the system, giving rise to the observed unidirectional motion. 
Second, we will carry out a detailed study of the ratchet dynamics of solitons governed by Eq. including not only 
the cases to = 2, our main focus in 28|, but also to = 4, which we compare to the previous one. Third, in order to 
show the generality of our results, we will present an analysis of the motion of the (/>^ kink, whose dynamics is given 

by 

0tt-0xx + 0^-0 = -/30t + /(t). (3) 

Finally, we will extend our previous work to analyze in depth the influence of dissipation in the directional motion. 
We will focus on the way that damping affects the phase values for which there is no soliton motion, and we will 
discuss the existence of values of the dissipation coefficient that lead to optimal soliton mobility. Interestingly, our 
results are reminiscent of the experimental findings in j2Cl |. where similar effects due to dissipation are found. 

The presentation of the above results is organized as follows: In the next Section we analyze twp nonlinear Klein- 
Gordon systems, specifically the sine-Gordon and (p'^ models, by means of a collective coordinate (CC) approach, 
obtaining the analytical and approximate expressions for the average of the velocity of the center of the soliton for 
TO = 2 and m — A. In Section III we verify our analytical results by means of numerical simulations. Section IV 
is devoted to the discussion of the most relevant effects of dissipation. To conclude the paper, in the last section 
we recapitulate the results of Sections III-V, make the connection with the experiments and summarize our main 
findings. 
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II. COLLECTIVE COORDINATE APPROACH 



A. Collective coordinate equations 



Our aim in this section is to obtain an approximate expression for the average vefocity of the center of the kink (as 
(f)'^ kinks are not sohtons, strictly speaking, we will refer to these excitations of both models as kink from now on). To 
this end, we will use one of the different CC approaches (see for a review). We introduce two variables, namely 
the center of the kink, X{t), and its width, l{t), through Rice's Ansatz [s^l, which amounts to specifying the kink 
solution as 



4>Q{x,t) — A: arctan 



exp 



x-X{t) 



for sG and 



(/)o(x, t) — tanh 



X ~ X{t) 



(4) 



(5) 



for 0^. Using the above expressions, and the variations of the energy and the momentum (two conserved quantities 
of the unperturbed systems Eqs. Q and Q) it can be shown (see j^ll for details) that the two collective variables 
obey the following two ordinary differential equations: 



P - 211- 2(311 = 



1 



1 

a 



where P[t) = MqIqX /l{t), fi/j = l/{^/alo) is Rice's frequency, and the parameters Mq, q, 
values according to the model, sG or (j)'^ (see Tabled). 



(6) 
(7) 

a and In take different 
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lo 
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7r7l2 


(tt" -6)/12 



TABLE I: Parameters corresponding to the two nonlinear Klein-Gordon models we work with, the sine-Gordon and 
q is the topological charge, Afo and lo are the mass and the width of the kink at rest. 



models. 



The quantity P{t) deserves some discussion on in its own. To begin with, it represents the momentum of the 
kink: Indeed, the expression for P{t) can be obtained as well by substituting Rice's Ansatz into the definition of 
the momentum P{t) = — dx (j)t(j>x- The same equations can be obtained using a projection technique with a 
Generalized Traveling Wave Ansatz (GTWA) (see details in H^)- In addition, it has been shown recently 33] that 
Eq. ((BJ for the momentum, obtained here within the framework of the Rice approximation, is an exact equation of 
motion, irrespective of the choice for the Ansatz. Hence, its solution is also exact and not an approximation. 

Let us now proceed with the analysis of Eqs. © and 0. After transients have elapsed and become negligible, 
[t':^ 1//3) the solution for Pit) is essentially given by 

P{t) = — Ve[ai sii\[5t + — Xi) + 02 sin(m(5t + 62— Xn 

where e = min(ei, 62) is a rescaling parameter, and 

J. 

Q 



Xi 



arctan 

q 



Xn 



arctan 



m5 

T 



and 02 



£2 



(8) 

(9) 
(10) 



We can now turn to Eq. (O, which is nothing but a parametrically driven oscillator: it is clear that l{t) is indirectly 
driven by the harmonic forces with frequencies 5 and m5, due to the term P^ on the r.h.s of Eq. Q. To solve 
(approximately) this equation, we expand l{t) in powers of e around the unperturbed kink width 



l{t) ^ lo + eh{t) + e^hit) + 



(11) 
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Substituting Eq. Hll() into Eq. ^ we find a hierarchy of equations for different orders of powers in e. The first three 
equations for li(t) are given by: 



hit) + ph{t) + nim = -^p'Wi + ^ + ^ (13) 

hit) + f3h{t) + nlut) = -ShP'm + ^ + ^ - Ifc - (14) 



?i(i)+/3/iW+l^^ji(t) = -^^'(*)'o, (12) 



These equations are hnear and therefore we can begin by solving the first one, Eq. (|12|l . by substituting the expression 
for the momentum, Eq. ©, into the r.h.s of (|12|l . With this, the equation for li becomes 

li{t) + (3h{t)nyi{t) = A1+A2 cos{2St + 201 - 2xi) + A3 cos{2mSt + 26*2 - 2xm) + 

+ cos[(m - l)5t + 02-91- (x™ - Xi)] - cos[(m + l)5t + 01+02- (xm + Xi)]}, (15) 

where 



Ai = -A2-A3, (16) 



- (17) 



-3 - i^ff' ^^^^ 

It is important to notice on the r.h.s of Eq. H15(l the presence of harmonics with frequencies 2S, 2mS and (m ± 1)6. 
As a consequence, all these frequencies appear in the expression for li(t) and they remain after a transient time 
t ^ 1//3, i.e., asymptotically we have 

, Ai A2sm{25t + 20i-2xi+02) A3sm{2mSt + 2e2-2xm + 02m) 



^(f^^- 452)2 +4/3252 ^(f]2j - 4m2,52)2 + 4m2/3252 



^4 sin[(m - l)5f + 6>2 - 6>i - (xm - Xi) + 

y/{n'j^ - (m - 1)2^2)2 + ^2(^ _ 1)2J2 

^4 sin[(TO + IjSt + 01+02- (Xm + Xl) + gm+l] 
^(172, - (to + 1)2^2)2 + ^2(^ 1)2^2 



(20) 



where 



0™ = arctan 



m(3S 



From the r.h.s of Eqs. Ijl^^ll and (|14|1 . the same reasoning as above leads one to find the harmonics corresponding to 
the second and third order corrections I2 (t) and I3 (t) . We do not include the explicit expressions for these corrections 
for the sake of brevity, but in Table^lwe collect all these values for to = 2, 3, 4. 



B. Velocity of the kink center 



Having found the exact expression for P{t) and the first terms in the perturbative expansion for l{t), we can now 
calculate the average velocity over one period T = 2Tr/S. To this end, we recall the previously defined expression for 
the momentum, P{t) — MqIqX /l{t). After transients have elapsed, which means after a time essentially given by the 
inverse of the dissipation coefhcient, the mean velocity of the kink can be expressed in terms of the CC as 



(21) 



5 



2"" harmonic((5) 


h 


h 


m 


2S, 2m5, {m±l)S 


25, 45, 4m5, (m ± 1)5, 
2(m ± 1)5, (m ± 3)5, (3m ± 1)5 


2 


S, 2S, 35, 45 


5, 25, 35, 45, 55, 65, 75, 85 


3 


25, 45, 65 


25, 45, 65, 85, 105, 125 


4 


25, 35, 55, 85 


5, 25, 35, 45, 55, 75, 
95, 105, 115, 135, 165 



TABLE II: Harmonic content of the first contributions to the perturbative expansion of l{t). Notice that 5 and m5 are the 
frequencies of the ac force (they appear in the momentum equation as well). 



Taking into account the expansion (CH), this expression can be written as 

^ . 1 rpimo+Mt)+e%it)+...) 

= (Xo(t))+e(i:i(t))+e2(i:2(t))+... (22) 

Therefore, as already stated, the mean velocity can be analytically (and approximately) calculated taking into account 
the exact solution for the momentum ||SJ| and the first terms of the expansion for the width of the kink. At order 
0{^), (Aro(t)) ^ {P{t)) = since the average of the momentum is zero [see Eq. ©]. Accordingly, the net motion 
of the kink can only arise from the contribution of higher order terms; hence, we proceed by computing the integral 
defining e{Xi{t)). Inserting Eqs. (jS)) and (^01) in (^3 . straightforward calculations show for m = 2 that 

. q^nlele2 ( 2 cos[20i - 9^ + (x2 - 2xi) - h] 



8M3 (/32 + ^2 ) ^p2^ 4^2 \^ ^(r!2j-<52)2+^2J2 

cos[2gi-g2 + (x2-2xi) + ^2] \ 

V(f72, - 4^2)2 + 4/32^2 I ■ 



(23) 



It is important to recall that for the CC approach to be valid, we must require that the kink shape is not very much 
distorted, and hence that ^ 1 and (5^1. Moreover, we would like to stress that expression H23|l is valid only for 
any finite but nonzero value of (3, since we have obtained it for t >> 

We can recast Eq. ()23|l in a simpler form to make more transparent the dependencies on the different parameters, 
as follows: 

e{Xi) = e2e2ri(/3,5)sin[02-20i + $i(/?,,5)], (24) 

where we have introduced an amplitude, Fi, whose specific form is not needed and is therefore omitted for brevity, 
and a phase. 



$i(/3, 6) = 2x1 - X2 - arctan 



2{nl - 5^)[{nl - 4(52)2 ^ 4^2j2] ^ (^2^ _ 4S^)[{ni _ ^2^2 + ^2^2] 



(25) 



which are only functions of the dissipation coefficient and the frequencies of the ac force. We note here that a similar 
expression, albeit starting from a less general Ansatz, has been obtained in |23| |. In this new form, it is apparent that 
the averaged velocity is a sinusoidal function of 0i and 02, and it is proporcional to e^e2. It is interesting to note 
that the same dependence on amplitudes and phases is also found in other systems, in principle unrelated with the 
ratchet dynamics of solitons we are discussing here (see j2l] | and references therein). On the other hand, from the 
experimental viewpoint, it is worth noting that, in the limit S <ti P, we find (V) ^ sin(20i — 62), in agreement with 
the experiments on J J reported in [T^. 

Proceeding with the analysis of the average velocity of the kink, from Eq. ()24|l it is obvious that net motion (i.e., 
with nonzero average velocity) is only possible if 

26li-e'2-$i 7^n7r,n==0,l,... (26) 

In this form, this condition, as will be shown below (and was pointed out in pol|'). arises from the breaking of the 
time-reversal symmetry by the biharmonic force. Indeed, $1 is related with the difference of phases between the 
harmonics with frequencies d and md in the functions P{t) and li{t). Thus, our first resonance criterion reads that 
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FIG. 1: a) Effective driving force in the sG model as a function of time, b) Trajectory of the center of the soliton; dashed hne: 
CC result obtained by numerically integrating Eqs. Q and 0; solid line: simulations of the full sG equation. Parameters used 
in both figures are S = 0.1, m = 2, = 0.05, 62 = -2.5, di = 0, ei = £2 = 0.2. 




time frequency/(2ji) 



FIG. 2; Ratchet dynamics in the sG model for m = 2. Left panel: Width of the kink vs time. Right panel: Discrete Fourier 
Transform (DFT) of the width of the kink. In both panels the dashed lines represent the numerical solution of Eqs. 16171 and 
the solid lines the computed width from the numerical simulation of Eq. (Q. Parameters are ei = e2 = 0.2, [3 — 0.05, S = 0.1, 
01 = -2.5 02 = 7v/2 - 2.5. 



the ratchet effect appears when l{t) oscillates with at least one of the frequencies of P{t); in addition, the difference 
of phases between the harmonics must be appropriate for the appeareance of net motion. In others words, fixing all 
the parameters of the ac force and damping except one of the two phases of the ac force, one can find a critical value 
of that phase that suppresses directed motion of kinks. The values at which the velocity vanishes are important, 
since for a force with a given frequency and damping, they indicate the phases at which the velocity of the kink is 
reversed. On the other hand, the changes in the phases are related with the breaking of the time-reversal symmetry 
of the system ^3^. 

Let us now turn to the case m = 4, a value for which the appearance of the ratchet effect is also predicted [23l |. 
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FIG. 3: Discrete Fourier Transform of the kink width. Left panel: m = 3 (no net motion); Right panel: m — 4 (ratchet motion 
of the kink takes place). Solid line: amplitude measured in simulations. Dashed line: amplitude as obtained by numerical 
integration of the CC equations ©-0- The parameters of the ac force and the dissipation coefficient are the same as in Fig. 

H 



In analogous manner to the preceding calculations, one can show that e(-'^i) is zero. Therefore, the frequencies in 
l{t) that contribute to the net motion must appear in higher order corrections (see Table llT)l . namely e^(X2)- After 
cumbersome calculations we find 



32Mo^(/32 + (52)2y(/?2 + 1652) 

/ sin(46li-6l2+X4-4xi+6'2-6ii) 



^/{nl - ,52)2 [3252^(^2^ _ 4^2)2 + 4^252 

(6(52 ^ cosj^ei - 02 + X4 - 4X1 - ^3 + 02-01) 



y^{n'j^ - .52)2 + f32p^^^2^ _ 4^2)2 + 4/3252^(^^2^ _ 9^2)2 + 9^2J2 

sin(46>i - 6>2 + X4 - 4xi - ^3 - ^1) 

^{nl - .52)2 + p2p^^Q2^ _ Qpy + 9^252 

sin(46'i - 6*2 + X4 - 4xi + O2 + h) 



2y^{nl - 452)2 + 4/3252^(^72^ _ 1652)2 + I6/3252 

(452 - r>2^) cos(4gi - 02 + X4 - 4x1 + 202 + 9i) \ 

' ^\(Sl\ - 452)2 + 4^2 J2] ^(f]2^ _ 1652)2 + 16/^252 J ' 



(27) 



which can be rewritten as 



^(X2) = e^2r2(/3,<5)sin[02-40i+$2(/3,<5)], (28) 



where again the amplitude, r2, and the phase, $2, are functions of the dissipation and the frequencies of the ac force. 
As for the mobility condition, which for the m = 2 case was given by Eq. H26() . a similar expression can be derived 
for the case m = 4 by imposing that the sine function vanishes in Eq. H28|l . 

For the case m = 3, the shift-symmetry is not broken and therefore ratchet motion can not take place. Indeed, the 
calculation of the average velocity based on the CC approach gives zero for all orders of the expansion. For this case 
the frequencies of the ac force (or the momentum) are odd harmonics of b {S and 3(5), whereas only even harmonics 
of 5 are found in the kink width oscillations (2n5, n E N). The complete selection rules for m = 2,3,4 appear in 
Table m where for a given m one identifies the harmonics corresponding to the oscillations of the width of the kink. 
Resonances take then place when at least one of these harmonics coincides with one of the two harmonics of the ac 
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FIG. 4: Dependence of the average velocity of a sG kink on the phase, 62, for different values of m. m = 2 (blue circles); m — 3 
(red diamonds) and m = 4 (black squares). Parameters are ei — €2 — 0.2, (3 = 0.05, S = 0.1, Oi = 0. 




0.1 0.2 0.3 



FIG. 5: Dependence of the function <I>i = —62 on the dissipation coefficient for m — 2 and /3 > 0.01. In the simulations as 
well as in the CC, we find 62 as the value of the phase between —n and for which the motion is suppressed. Solid line: CC 
result; dashed line: Eq. H25^ : symbols: numerical simulations of the full sine-Gordon equation. Parameters are ei = £2 = 0.2 
and 6*1 = 0. 

force. We can verify that if this condition is fullfiled, the shift symmetry is broken and viceversa. In the simplest 
case, m = 1, for which the driving consists only of a single harmonic, the ratchet effect is also absent. Looking at 
Eq. Q it is immediately realized that only terms with frequency 26 appear in the equation Q , so that the resonance 
condition is not achieved in agreement with previous results '31i|. This discussion in terms of resonant frequencies 
can be extended, in an analogous but cumbersome calculation, to any positive integer number of the frequency of the 
second harmonic, i.e., for higher integer values of m. Finally, an interesting remark is that the resonance condition 
can also be interpreted as a synchronization between the ac force and the oscillations of the width of the kink, i.e., 
directional transport occurs whenever the ac force is locked to the oscillations of the width ■ 
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C. Mechanism for directional motion 



In the previous subsection we have estabhshed the conditions for net motion of the kink to take place. We argued 
that the fact that the width of the kink is a dynamical variable is not a sufficient condition for the existence of directed 
motion and, on the contrary, net motion can arise only when, first, at least one of the two harmonics of the ac driven 
force appears in the oscillations of the width of the kink (which is equivalent to breaking the shift symmetry) and, 
second, when the time-reversal symmetry is broken. However, these two conditions, albeit simply stated, do not allow 
for a physical or intuitive understanding of the mechanisms at work in the system. In order to shed light on those 
mechanisms, it is convenient to write Eq. ^ as 

Cast in this form, it is easy to realize the existence of an effective driving force, given by 

whose amplitude is modulated by the oscillations of the width. The consequence is that, while the average of the 
external ac force is always zero, the average of this effective force term vanishes when the harmonics of the external 
driving force f{t) do not match those in the oscillations of the width or when they oscillate with at least one common 
frequency but with appropriate critical phases. Such a nonzero effective driving force is in fact the same as having a 
dc force acting on the kink. This means that the energy is inhomogeneously pumped into the system |3(i| |. through a 
process that is modulated by the kink width oscillations. 

On the other hand, Eq. H29|l has a term XI /I that describes the energy transfer from the traslational mode to the 
internal mode. Instead of analyzing which portion of energy delivered by the ac force is used for the dynamics of the 
internal mode, we find it more illuminating to study how different this transfer is in the two halves of one period. 
To this end, Fig. da) presents the effective force for the case m = 2 for certain parameters. The difference between 
the brown and grey areas divided by the period represents the average of the effective force in one period. The 
computation shows that (fef) — 0.0055. Consequently, the center of the kink should move in the positive direction 
and, indeed, as is depicted in Fig. ^b), such motion takes place. In addition, the mean velocity is (v) = 0.108, 
the kink moves with constant velocity on average, and hence the average acceleration is zero. One can check that 
(3 (v) — 0.0054 in close agreement with the result found above for the effective force. We also computed the average 
of the coupling term which gives {XI /I) — —0.0001. This tells us that, on average, the coupling term is equivalent to 
a friction force but also that its value is much smaller compared to the average of the driving force. Hence, most of 
the net force is used in the traslational motion. 



III. NUMERICAL SIMULATION RESULTS 



The analytical results of the previous section were derived within the CC approach, which amounts to consider 
that perturbations affect only the X and I variables of the kink but not its shape, i.e., phonons have been neglected. 
Therefore, we expect that the results will be valid for small amplitudes and frequencies of the ac force, and also not too 
small damping coefficient. In any event, they must be compared to numerical simulations of the full partial differential 
equations JQ) and To this end, we have integrated numerically both equations, using the Strauss Vazquez scheme 
[33 as well as a fourth-order Runge-Kutta method, choosing a total length of L = 100,300, with steps At = 0.01, 
Aa:; = 0.1. We have used aperiodic boundary conditions with a kink at rest as initial condition. In the following we 
will monitor the dynamics induced by the ratchet effect through the dependence of the average velocity on the phases 
of the bi- harmonic force and the dissipation coefficient, which in turn requires the computation of the center and the 
width of the kink from the simulations of the full equations. 

For the numerical calculation of the center and width of the kink we have improved the procedure suggested in |38j | , 
taking into account the oscillations of the ground state due to the action of the ac driving. In particular, for the sG 
kink (and for the (p'^ kink, changing parameters as indicated), this method reduces to searching, for a given and fixed 
time, in the discrete lattice, those points a;„ and Xn+i such that (pn ^ + (t>va,c{t) and (j)n+i > n + 0vac(*), where 
(pn = 4>{xn,t), (pvacit) = ~^ ±00, t) represents the oscillations of the background field. In this case this function 
can be computed as 4>n — 27r (0jv — 1 for 4''^) where = L/ Ax, L being the total length of the numerically simulated 
system. Subsequently we estimate, by using a linear interpolation method with two points {xn, 4>n) and {xn+i, 4>n+i), 
the corresponding point Xn (it will be our computed center of the kink X{t)) where = tt + (/ivac {4> — 4'vac for 
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FIG. 6: Dependence of the average velocity on 62 for m = 2. Left panel: Parameters are ei — e2 — 0.02. Right panel: 
ei — 62 — 0.2. Other parameters are P — 0.05, S — 0.1. In both panels two values for 9i are considered: blue squares are 
the results of the numerical simulations for 9i — 0, and black circles correspond to 9i = 7r/2. Solid lines correspond to the 
numerical solution of the CC Eqs. JHJ-Q. In the left (right) panel dashed lines correspond to {e{Xi) /5), obtained from 

Eq. 



Second, in order to compute the width of the kink, we look for the value of l(t) that minimizes the expression 



k=l 











<pk - 1 


00 


[ m \ 


+ 0vac^ 



where 0o is defined by using Q (or (O in the case of the (f)^ system) and fc = 1, 2...n...N. 

In the following subsections, we will present the numerical analysis in two parts according to the two models we 
deal with, namely sG and (p'^. In each system we will analyze the ratchet dynamics of the kink for 771 = 2,4 and the 
zero average velocity for m — 3. 



A. Sine-Gordon model 



Our first comparison between the CC analyisis and the numerical simulations is presented in Fig. ^b), where 
the time evolution of the center of the kink is presented. As can be seen from the plot, the agreement between 
our analytical approximation and the numerical results is excellent, even for long times, suggesting that we can be 
confident that the CC approach captures the ratchet dynamics of the system. Nevertheless, to complete the analysis 
and fully confirm our analytical predictions, it is crucial to understand how the width of the kink evolves when it 
is driven by the bi-harmonic external force. The evolution in time of the width of the kink as obtained from the 
numerical simulations of the full sG equation and from the CC equations are compared in the left panel of Fig. |21 
Once again, we observe an excellent agreement between simulations and the CC framework. In order to validate 
our predictions about the resonance criterion, we proceed to the determination of the Fourier components for the 
oscillations of the width of the kink (see right panel of Fig. The Discrete Fourier Transform (DFT) shows an 
excelent agreement between the numerical simulations of Eq. and the numerical solutions of Eqs. (|6I7|) . confirming 
the resonance criterion. For instance, we can observe that, for m — 2, the spectrum of lit) contains the frequencies 
5 and 25, the same frequencies of the bi-harmonic external force. Consequently, one could expect a net motion of 
the soliton, as was pointed out in Table ^ and as the numerical simulations show. To further confirm the resonance 
criterion, we analyze the spectrum of l{t) for m — 3,4, plotted in Fig. |31 For m = 3 neither 5 nor 3(5 appear in the 
spectrum, thus confirming our predictions (see Table ITl|: We observe only the presence of even harmonics, and hence 
we expect an oscillatory motion of the center of the kink similar to that obtained for a force with only one harmonic 
component |26|. Conversely, for m = 4 we observe the appearance of frequencies 5 and AS. Therefore, the occurrence 
of a net motion is also expected as in case of ?ti = 2 and is indeed confirmed by the simulations (not shown). 



11 




-0.03 



FIG. 7: Dependence of the average velocity on 02 for m = 4. Circles correspond to numerical simulations; solid line to the 
numerical solution of CC Eqs. ©-Qi and dashed line stands for e^{X2)/9. Parameters are ei = e2 — 0.2, 9i = 0, S — 0.1 and 
P = 0.05. 



Having verified our resonance criterion for the ocurrence of ratchet phenomena, we now turn to discussing how the 
average velocity depends on the parameters of the ac force and damping. In Fig. ^ we plot the average velocity of the 
kink, computed from the numerical simulations of the full PDE, for m = 2, 3, and 4 versus the phase of the second 
harmonic 02 [see Eq. (j^J]. In this plot, the predictions for the existence or not of motion for different values of m are 
confirmed. Notice the sinusoidal dependence of the mean velocity on the phase, as predicted by the expressions H23|) 
and H27() . Another relevant feature of Fig.^jis the difference of the mean velocity between the cases m = 2 and m — A. 
In principle, net motion can occur for any even value of m; however, we observe that the maximum of the velocity 
decreases when m is increased. This can be understood from Eqs. (|24|l and (|28|l . which predict that {V) ^ e\e2 for 
m = 2 and {V) ~ e\e2 for m — A. For greater values of m, {V) ^ e™e2 21], so if ei and £2 are small enough, the 
velocity goes to zero as m is increased. Furthermore, the amplitudes of the peaks of the DFTs corresponding to the 
frequencies 5 and 25 for m = 2 are greater than their counterparts for m = A, d and AS (see right panel of Figs. [21 
and|21l. In this figure we also notice that there are some values of the phase 62 for which there is no net motion: For 
example, for m = 2, these values are approximately given by 62 — —0.9 + mr (n = 0, 1, ...). We will refer to those 
values of the phase for which there is no net motion as critical values. For m = 2, the whole set of critical phase values 
as a function of dissipation and frequencies is provided by the expression 62 = — <I>i(/3, S) + mr. These correspond to 
those values of the phase for which the condition H26II is not fulfilled. 

The dependence of 6) on the dissipation coefficient and the frequency of the ac force is shown in Fig. The 

figure shows that for the largest value of the frequency, the agreeement between the CC theory and the simulations of 
the full partial differential equation is only qualitative. The reason for this lies in the fact that, when the harmonics of 
the ac force contain the frequencies S and 2S, the width of the kink is excited with frequencies 26 and AS, respectively 
(because of the term in see also previous studies ■ail ill). If we chose S = 0.25, l{t) wiU osciUate with 
2S = 0.5, AS = 1, or even higher frequencies, which are inside the phonon band, which lies above unit frequency. 
Therefore, for large enough amplitude of the force and small dissipation, strong excitation of phonons is expected 
(even kink-antinkinks can appear) and correspondingly the failure of our CC approach, since it does not take into 
account the phonon contribution to the motion. This is yet another reason why simulations for very low values of the 
dissipation coefficient (and hence with larger phonon production) are not considered in the plot. 

Let us now concentrate on the m = 2 case and compare the average velocities as a function of 62 , computed from 
the direct numerical simulations of Eq. ^ and from the numerical solution of the CC equations ©"tZI)- In Fig- 
|B1 we observe an excellent agreement between these two average velocities. In the left panel, overimposed to these 
values, we show also the approximate values of (V) obtained from Eq. I|23|) . In the right panel, the amplitudes of 
the two harmonics have been increased by one order of magnitude; this leads to noticeable phonon production in 
the simulations and consequently to a quantitative disagreement between the simulations and the CC approach, and 
therefore a factor 1/5 has been introduced in order to adjust the values of {V) obtained from Eq. (|23|) to the results 
of the numerical simulations and the numerical solution of CC equations (as is obvious, the qualitative agreement is 
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FIG. 8: Average velocity of a sG kink as a function of the dissipation for P > 0.02. Panel (a): ei — e2 — 0.02, 5 — 0.1, 9i — 7r/2; 
Numerical simulations: 62 ^ n {*); 02 = (+). Panel (b): ei = 0.04, £2 = 0.026, S = 0.25, 61 = -7r/2, 82 = -tv/2 + 0.8; 
Numerical simulations: (□). Panel (c): ei = 0.04, £2 = 0.026, 5 = 0.25, 9i — —n/2, 62 = — 7r/2; Numerical simulations: {()). In 
each figure the solid line shows (V) obtained from the numerical solution of Eqs. 0-0, and the dashed line is {V) calculated 
from Eq. CT . 



perfect and in particular for the critical phase values the agreement is quantitative). For the case m = 4 the situation 
is the same: There is good agreement between the results of the CC equations and the simulations, although once 
again the analytical expression {V) fits only qualitatively to the numerical results (see Fig.O. 

One important point, relevant in experimental contexts, is the prediction of the CC theory about the dependence of 
the mean velocity on the dissipation coefficient. According to the standard behavior of point particles under friction, 
one would expect a monotonic decreasing of the average velocity to zero as a function of the damping. However, 
the existence of an optimal damping for the occurrence of net motion, and the appeareance of current reversal upon 
varying the damping (see O), show a richer and more complex behaviour in the case of the soliton ratchets 
we are discussing, depending on the parameters of the ac force. These phenomena, expected also from Eqs. (|23|l and 
(|27|l . are shown in Fig. |S1 where for the first sets of parameters (Fig. we observe that the velocity drastically 
decreases to zero as (3 is increased. Notice that in this case the time-reversal symmetry of the ac force is broken in 
the optimal way [2]| by the given phases. In Fig. IS]:, when the phases are chosen in a manner such that the time 
reversal symmetry f{t) = f{—t) is fulfilled in the absence of dissipation, an optimal value of /3 for the transport is 
observed. In this case, when the dissipation coefficient increases starting from (3 = 0, the time-reversal symmetry 
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FIG. 9: Average velocity of a sG kink as a function of the dissipation for 0.01 < l3 < 0.45 and for two different frequencies: 
S = 0.1 (squares); S = 0.25 (circles). The other parameters are ei = 0.2, £2 = 0.12, 9i — —n/2, 02 = — 7r/2. In both cases the 
solid lines show the results obtained from the numerical solution of Eqs. @-(I7J. 

begins to break, and therefore an increase of the average velocity is expected. On the other hand, increasing the 
damping tends to suppress any motion. Hence, in between these two opposite mechanisms an optimal value of the 
damping must appear ,21.] . Between these two limit cases, we can see the appeareance of current reversal (see Fig. 
|H1d) as was pointed out in p^ . 

From Fig. |Slwe can see that the best agreement between CC theory and numerical simulations is obtained for small 
values of the frequencies of the ac force. Indeed, in Fig. for 6 = 0.25, due to the phonon contribution (see the 
discussion above), we observe only a qualitative agreement between the numerical simulations and the CC theory. 
This disagremeent is even clearer in Fig. where we take large values of the amplitudes of the ac force and compare 
the mobility of the kink for 6 — 0.1 and S — 0.25. We can observe that for the largest value of frequency together with 
the smallest values of /3 the results from the theory do not fit quantitatively the values obtained from the numerical 
simulations. Furthermore, we even observe current reversal in the results of the numerical simulations, contrary to the 
predictions of the CC theory. As we explained above, this disagreement arises from the large production of phonons 
for this value of the frequency, unaccounted for in the CC approach. 

B. model 

After describing in the previous subsection the ratchet dynamics of kinks in the sG equation, here we focus our 
attention on the (p* model. This is important in order to ascertain to what extent does the phenomenology found 
depend on the unperturbed equation being integrable, or on the structure of the sG soliton (which, for instance, 
does not have internal modes |39j, whereas the (j)'^ model does have one). In view of this, we carried out numerical 
simulations of the (j>'^ model, to compare with the analytical CC results in Sec. ^ As we will now summarize, the 
phenomenology of directed motion of kinks in the 0"* model is the same as in the sG model. 

To begin with, we can observe in Fig. IIUI that our first resonance criterion on the synchronization between the 
oscillations of l{t) and f{t) does also apply for the model. As we can see, for even value of to (to = 2,4) there are 
resonant frequencies, while in case of m = 3 the frequencies of the ac force S and 3 S never appear in the spectrum of 
the width of the kink. Correspondingly, in our numerical simulations we have found unidirectional motion only in cases 
of m = 2, m = 4, in contrast to to = 3 where an oscillatory motion of the center of the kink takes place (not shown). 
The sinusoidal dependence of the average velocity on the phases is found as well as in the sG equation. The other 
feature we focused on, namely dissipation effects, agrees with the previous picture as well: In Fig. llll we observe, as in 
the sG case, that changing the damping we can decrease the average of the velocity, rectify the motion, or optimize it. 
In this respect, another interesting issue is the comparison between the observed mobility of the kink between the sG 
and (j)'^ models. From Eq. H29I) . it can be noticed that, if we consider the first approximation of l{t), the center of the 
kink will feel an ac force with the amplitudes of the two harmonics modulated by the factor A ~ q/Mq, which in sG is 
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FIG. 10; Discrete Fourier Transform of the width of a (f)'^ kink, (a): m = 2; (b): m = 3; (c): m = 4. Solid line: Simulations of 
Eq. @. Dashed line: numerical integration of the CC equations. 

XsG = 0.785 and in (f)^, = 2.121. Therefore, for the same values of e.^ in both systems, the effective amphtude of 
the ac force that the center of the kink experiences in (j)'^ is greater than in sG, Xsc^f^ < X^^ef'^. Therefore, from the 
CC analysis we expect higher mobility for the 0"* model, which is in good agreement with the numerical simulations 
as is shown in Figs.lSlandllll We thus conclude that the existence of net motion of kinks (topological excitations) is a 
generic phenomenon, both in terms of the criterion for its appearance and its main characteristics, and not a specific 
feature of the sG equation. 

IV. CONCLUSIONS 

In this paper we have investigated in detail the behavior of topological excitations when driven by ac forces in 
homogeneous systems. We have focused on the nonlinear Klein-Gordon family of models as a paradigmatic example, 
and specifically on the sG and cj)'^ equations as examples of integrable and non-integrable models, respectively. By 
means of an analytical technique of the CC class, we have shown that two conditions must be fulfilled in order to 
observe unidirectional motion of the kinks: First, the shift symmetry must be broken, and second the time-reversal 
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FIG. 11: Average velocity of a (/!>^ kink as a function of the dissipation for (3 > 0.02. (a): ei = £2 = 0.02, S — 0.1, 6i = it/2. 
Numerical simulations: 02 = tv (V); 62 = (a). (b): ei = 0.04, £2 = 0.026, d = 0.25, Oi = -n/2, 62 = -it/2 + 0.8. Numerical 
simulations: (+)• (c): ei = 0.04, £2 = 0.026, 5 = 0.25, 9\ = —n/2, 62 = — 7r/2; Numerical simulations: (x). In all cases the 
solid lines show results obtained from the numerical solution of Eqs. JSJ-Q and the dashed lines are the values obtained from 

Eq. (unj. 



symmetry must be broken as well. In the absence of dissipation, the first symmetry may be violated by a bi-harmonic 
force containing even and odd harmonics of the frequency. We have devoted a large part of our work to characterize 
the case in which the force is composed of the first and second harmonics [m = 2 in Eq. m], and we have found 
that the velocity of the net motion depends on the relative phase between the two harmonic forces in a sinusoidal 
manner. The values of the phase for which the velocity vanishes (critical phase values) are precisely those for which, 
loosely speaking, the breaking of the time-reversal symmetry by the ac force is compensated by the breaking of the 
time-reversal symmetry due to the damping. Numerical simulations for the two models considered are in very good 
agreement with the analytical results, qualitatively in terms of the velocity value, and quantitatively as far as the 
critical phase values are concerned. We have also studied the cases of the third (m = 3) and fourth (to = 4) harmonics, 
finding, respectively, no motion (because the shift symmetry is not broken) and motion slower than in the previous 
case. Indeed, as we have seen, in general the efficiency of the driving diminishes when the second harmonic force is a 
higher harmonic. In all cases, the arguments in terms of shift symmetry breaking forces can be equivalently presented 
as existence or not of resonant behavior between the frequencies present in the driving and those in the time evolution 
of the width of the kink. 

In addition to the breaking of the time-reversal symmetry of the ac force, we have also studied another form of 
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time-reversal symmetry violation, namely by introducing dissipation in the models. In the presence of dissipation, 
we have observed three different types of behavior depending on the parameters of the driving: A monotonically 
decreasing velocity of the kink with increasing dissipation, the appearance of an optimum value of the dissipation 
with a maximum in the velocity, and even reversals of the directed motion for some values of the dissipation. In 
all three cases, the CC theory describes correctly the phenomena as observed in the simulation, except when the 
dissipation is very small, allowing for a non-negligible excitation of phonons in the system that breaks down the basic 
assumption of the CC theory. 

The conclusions we have just summarized acquire greater relevance when compared to different types of experimental 
results [TEl ITM l20j . In the absence of dissipation , experiments in optical lattices with biharmonic driving find a 
sinusoidal dep endence on the relative phase similar to the one reported here. When dissipation is present in Josephson 
junctions [131, the same sinusoidal dependence is again found, and optimal driving parameters appear. Finally, recent 
work on dissipation effects on optical lattices exhibits several of the features we have just discussed, such as the 
dependence of the critical phase on the damping, and the overall sinusoidal dependence. Although the motion of 
cold atoms is essentially based on a one-particle description the dissipation in the system is caused by ground state 
transitions, which manifests the presence of an internal dynamics. The role of this internal dynamics for the ratchet 
mechanism in cold atoms has been put forward in 40] . In this regard the role of the internal mode in our CC approach 
could shed light on the effect of the interband transitions and their disipative effects in the ratchet transport of cold 
atoms. We thus sec that our results are very generic, in so far as they apply to a large family of theoretical models, 
both integrable and non-integrable, and they capture the essential ingredients of experiments in different fields 2^. 
Therefore, the main conclusion of our work is that we have correctly identified the mechanism for the appearance 
of directed motion due to zero-mean forces through the breaking of symmetries, and we have provided a physical 
interpretation of this mechanism in terms of the coupling between excitation modes of the system. We hope that this 
work encourages further experimental work to check the remaining features discussed. 
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